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Abstract: Numerical simulation of lattice QCD is now an indispensable tool for investi¬ 
gating non-perturbative properties of quarks at low energies. In this paper, I introduce a 
recent lattice QCD calculation of light hadron spectrum and quark masses performed on 
a dedicated parallel computer, CP-PACS. From a simulation of QCD with two flavors of 
dynamical quarks, we hnd mJ^®(2GeV) = 3.441 q; 22 MeV for the averaged mass of u and d 
quarks to reproduce the experimental masses of tt and p mesons, and m^®(2GeV) = 88tg 
MeV or QOln MeV to reproduce the A or 0 meson mass. These values are about 20-30% 
smaller than the previous estimates using the quenched approximation. 


1 Introduction 

Quarks are among the most fundamental building blocks of nature. Unlike the previous 
elementary constituents, quarks are found to possess a remarkable property, confinement, 
he., we cannot isolate quarks by dissociating hadrons. Historically, fundamental theories 
of nature have emerged through reactions of isolated elements to weak external perturba¬ 
tions. Quark confinement, however, has the consequence that we cannot directly measure 
the fundamental properties of quarks by experiment. 

Fortunately, quarks are found to possess another remarkable property, asymptotic free¬ 
dom, whereby interactions among quarks become weak at high energies. This enabled us 
to identify quantum chromodynamics (QGD) as the fundamental theory of quarks. QGD 
has been quite successful in explaining high energy processes of particles. 

On the other hand, because the coupling parameter of QGD becomes large at low en¬ 
ergies, basic properties of hadrons, such as mass spectrum, spatial sizes, decay constants, 
etc., are not calculable by conventional analytic methods. In other words, we cannot 
reconstruct hadrons from quarks based on the dynamics of QGD yet. In order to answer 
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the fundamental question of whether QCD is correct also at low energies, we have to carry 
out non-perturbative calculations of hadrons directly from the hrst principles of QCD. 

Calculation of low-energy hadron properties is also required for determination of fun¬ 
damental parameters of nature. Due to quark conhnement, quark masses and coupling 
parameters have to be inferred indirectly from a comparison of experimental results for 
hadron masses etc. with a theoretical calculation of them as functions of the fundamental 
parameters of QCD. Also for electroweak interactions, although the core part of the reac¬ 
tion can be reliably calculated in perturbation theory, uncertainties in QCD corrections 
to the quark currents propagate to the determination of the Cabibbo-Kobayashi-Maskawa 
parameters. 

Numerical simulation of quarks based on the lattice formulation of QCD is currently 
the only method to calculate low-energy non-perturbative properties of hadrons. Re¬ 
cently, predictions with big impact are beginning to be produced through development of 
dedicated parallel computers |l[]. In particular, a big breakthrough was achieved by the 
CP-PACS computer developed at the University of Tsukuba 0. In this paper, I summa¬ 
rize major results from CP-PACS, focusing on the most fundamental topics of QCD, he., 
the light hadron spectrum and light quark masses. 

In Sec. 1^, I introduce the lattice formulation of QCD and explain why dedicated com¬ 
puters have been developed by physicists. Results for the light hadron spectrum and light 
quark mass are presented in Secs. and Conclusions are given in Sec. 

2 Lattice QCD and dedicated parallel computers 

We formulate QCD on a four-dimensional hypercubic lattice with a finite lattice spacing a 
and finite lattice size Na. The real world is defined by the limit of vanishing lattice spacing 
a ^ 0 keeping the lattice size Na finite or sufficiently large (the continuum limit). Before 
taking this limit, the theory is finite and mathematically well-defined. Therefore, we can 
apply various non-perturbative techniques. Here, we perform numerical simulations on 
finite lattices to calculate hadrons. 

To accomplish this, the most time-consuming part of the calculation is the inversion 
of quark propagation kernels. Numerically, this is an inversion of large sparse complex 
matrices. For the case of Wilson-type lattice quark actions, the typical size of the matrix 
is 12U X 12U, where 12 = 3 x 4 is the freedom of color and spin, and V = N‘^ is the lattice 
volume. Because the condition number is inversely proportional to the quark mass, the 
inversion is numerically more intensive when we decrease the quark mass. Even with the 
latest supercomputers, it is difficult to simulate the light u, d quarks directly. Therefore, 
in addition to the continuum extrapolation, we need to extrapolate the results at typically 
around the s quark mass to the physical u, d quark mass point. This procedure is called 
“the chiral extrapolation”. 

In summary, in order to extract a prediction for the real world, we have to perform 
continuum and chiral extrapolations. To get a precise and reliable result, it is essential to 
have good control of these extrapolations. This requires a large-scale systematic calcula¬ 
tion, and thus huge computer power. The power of vector supercomputers in 1980’s was 
not sufficient, which motivated several groups of lattice physicists to construct parallel 
computers dedicated to lattice calculations |ll]. 

The historical development of computer speed is shown in Fig. |l|. The frontier of 
speed has been advanced by parallel computers since around 1990. Dedicated machines 
developed by physicists contributed much to this trend. At the University of Tsukuba, we 
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Figure 1: Development of the computer speed. 


Figure 2: Quenched light hadron spectrum |^. 


have developed two machines, QCDPAX (1990) 0 and CP-PACS (1996) 0. CP-PACS 
is a parallel computer achieving peak performance of 614.4 GFLOPS with 2048 single¬ 
processor nodes. Since the hrst power-on in 1996, intensive lattice QCD calculations have 
been made on CP-PACS. In the following, I focus on the studies of hadron spectrum and 
light quark masses. 


3 Light hadron mass spectrum 

Precise calculation of the hadron mass spectrum directly from the first principles of QCD 
is one of the main goals of lattice QCD. Because the computer power required is enormous, 
we study the issue step by step as follows: 

Step 1 : Calculate in the quenched approximation, in which the effects of dynamical pair 
creation and annihilation of quarks are neglected. 

Step 2 : Include the dynamical u, d quarks in a degenerate approximation, while heavier 
quarks are treated in the quenched approximation {two-flavor full QCD). 

Step 3 : Include the dynamical s quark {2+1 flavor full QCD). 

Steps 4, 5, ■ • ■ : Introduce the u, d mass difference, dynamical c quark, etc. 

3.1 Quenched studies 

With the quenched approximation, we can reduce the computer time by a factor of sev¬ 
eral hundred preserving the basic properties of QCD (conhnement, asymptotic freedom, 
and spontaneous breakdown of chiral symmetry). The effects of the approximation are 
expected to be about 10% in the spectrum. 

The hrst studies in this approximation had already been made in early 80’s. The 
issue turned out to be quite tough and computationally demanding. Actually, it took 
about ten years until the hrst systematic study, performing all the extrapolations, was 
attempted in 1993 [Q. From this study, the quenched light hadron spectrum was found 
to be consistent with experiment within the errors of about 10%. However, the quality of 
data was insufficient to test the accuracy of numerical extrapolations, and also the hnal 
errors were too large to resolve the quenching artifact. 

The status was signihcantly improved by CP-PACS . From an intensive computation 
involving about 100 times more hoating-point calculations than the previous studies, it 
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Figure 3: Continuum extrapolation of vector meson 
masses and Mk* in two-flavor full QCD [Nf = 
2) and quenched QCD (TV/ = 0), using Mk as input 
0 . Open circles are recent results using a different 
improved lattice action [|j. 


Figure 4: Continuum extrapolation of the s quark 
mass TOs in the MS scheme at 2 GeV from two- 
flavor 0 and quenched QCD ||] . Quenched masses 
are shown by triangles. Filled and open symbols 
are for the Mk and inputs, respectively. 


became possible to perform systematic tests on the quality of the extrapolations. Errors 
in the hnal results for the light hadron spectrum are now conhdently estimated to be less 
than about 1% for mesons and about 3% for baryons. These errors include statistical and 
all systematic errors except for those from the quenched approximation itself. 

The quenched light hadron spectrum is shown in Fig. |^. Experimental values of meson 
masses — Mp, and either Mk or — are used as inputs to £x the lattice spacing a 
(or the strong coupling constant), the average u, d quark mass rriud, and the s quark mass 
rUs- Masses for other hadrons are among the predictions of QCD. From Fig. we hnd 
that the global pattern of the spectrum is well reproduced. At the same time, we clearly 
see that the quenched spectrum deviates from experiment by about 10% (7 standard 
deviations) for the worst case: Results from the Mjy-input and M^-input are discrepant 
by, and the hyperhne splitting between K and K* mesons is smaller than experiment by 
about 10%. Decuplet baryon splittings are also small. Because all other systematic errors 
are well controlled, we identify the discrepancy as quenching artifact. 

3.2 Two-flavor full QCD 

From the quenched simulation we And that, to calculate hadronic quantities with precision 
better than 10%, we have to incorporate dynamical quarks. A naive extension of the 
quenched simulation to full QCD is difficult because several hundred times more computer 
time is required. A partial solution is given by improvement of the lattice theory, whereby 
continuum properties are realized on coarser lattices. From a preparatory study 0, we find 
that the combination of renormalization-group improved glue action and clover-improved 
Wilson quark action is effective in removing major lattice artifacts. This reduces the 
computer time by a factor of about ten. 

Although more might still be hoped for, we can start the hrst systematic studies of two- 
flavor full QCD on the CP-PACS performing both chiral and continuum extrapolations 0. 
To identify dynamical quark effects clearly, we carried out another quenched simulation 
using the same improved action. 

Figure]^ shows the continuum extrapolation of vector meson masses. We And that the 









Table 1: Light quark masses in MeV in the MS scheme at 2 GeV [Q. 
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two quenched results (open and filled squares) lead to universal values in the continuum 
limit a = 0. They deviate from the experimental values (diamonds) as noted in the 
previous subsection. On the other hand, the full QCD results (circles) extrapolate to 
values much closer to experiment. Accordingly, we find no big discrepancies between the 
results from and inputs. This means that the quenching artifacts in the spectrum 
are mostly removed by introducing dynamical u, d quarks. 

The whole light hadron spectrum can now be approximately reproduced from QCD by 
adjusting just three parameters; the strong coupling constant, mud and mg. This provides 
us with a strong confirmation that QCD is correct also at low energies. Remaining small 
deviations from experiment may be explained by the quenched approximation of the s 
quark. In order to confirm this, however, uncertainties from the chiral and continuum 
extrapolations should be reduced to the level of our quenched study. This goal is reserved 
for future work. 


4 Light quark masses 


Adjustment of mud and in the calculation of the hadron spectrum provides us with the 
most direct determination of quark masses from QCD. Fig. ^ shows mg in the MS scheme 
at /i = 2 GeV as functions of the lattice spacing a. Results from the quenched improved 
action, which are consistent with the quenched standard action in the continuum limit, 
are omitted for clarity. 

On the lattice, there are several alternative definitions for the quark mass, with the 
difference being 0{a). In Fig. ^ they are denoted as VWI and AWI (vector and axial- 
vector Ward-Takahashi identity). See for details. As shown in Fig. different 

definitions lead to different values of mg at finite lattice spacings. This has been a big 
source of uncertainty in the previous calculations. From Fig. we see that they converge 
to a universal value in the continuum limit, as theoretically expected. This confirms the 
quality of our calculations. 

The quenched value for mg, however, differs by about 20% between Mi^-input and 
M<^-input. This is equivalent to the quenched artifact discussed in Sec. We note that 


the quenched artifact is larger in the light quark masses than in the hadron spectrum. 
When we turn on the dynamical u, d quarks, the discrepancy between the inputs disap¬ 
pears within our errors, in accord with the fact that the full QCD reproduces the hadron 
spectrum better (Sec. [^) . 

The light quark masses in the continuum limit are summarized in Table |^. The errors 
include our estimates for systematic errors from chiral and continuum extrapolations and 
renormalization factors. We note that the masses from two-flavor QCD are 20-30% smaller 
than those from the quenched QCD. In particular, our s quark mass in Nf = 2 QCD is 
about 90 MeV, which is significantly smaller than the value ~ 150 MeV often used in 
phenomenology. Our results are, however, consistent with recent estimates mg = 83-130 
MeV and mud = 3.4-5.3 MeV from QCD sum rules H and mg/mud = 24.4 ± 1.5 from 
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5 Conclusions 


An intensive calculation of the light hadron spectrum in the quenched approximation has 
revealed discrepancies of about 10% from experiment. With two flavors of dynamical u, d 
quarks, the hrst systematic study performing both the continuum and chiral extrapola¬ 
tions has shown that the discrepancies are mostly removed, providing us with a strong 
conhrmation that QCD is the correct theory of quarks at low energies as well. 

Precision calculation of hadron masses enables us to determine fundamental parameters 
of quarks directly from the hrst principles of QCD. We hnd that the dynamical quark 
effect is as large as 20-30% in light quark masses. Noticeable and sizable dynamical quark 
effects are observed also in B meson decay constants pd|, the equation of state at high 
0 : 


and the topological structure of the QCD vacuum |]T3. At the same 


temperatures 

time, new types of lattice fermions have begun to be applied to study hadronic matrix 
elements relevant to the AJ = 1/2 rule and the CP violation parameters |T^ . 

Because a shift in fundamental parameters can have signihcant implications to phe¬ 
nomenological studies of the standard model, it is urgent to evaluate dynamical quark 
effects in other hadronic quantities as well. The inhuence of dynamical s quark should 
also be studied. Studies in these directions are under way to clarify the precise structure 
of the standard model. 

The studies presented in this report have been performed by the CP-PACS Collabora¬ 
tion. I thank the members of the Collaboration for discussions and comments. This work 
is in part supported by the Grants-in-Aid of Ministry of Education, Science and Culture 
(Nos. 12304011 and 13640260) and JSPS Research for Future Program. 
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